Transmission properties in waveguides: An optical streamline analysis 
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A novel approach to study transmission through waveguides in terms of optical streamlines is 
presented. This theoretical framework combines the computational performance of beam propa- 
gation methods with the possibility to monitor the passage of light through the guiding medium 
by means of these sampler paths. In this way, not only the optical flow along the waveguide can 
be followed in detail, but also a fair estimate of the transmitted light (intensity) can be accounted 
for by counting streamline arrivals with starting points statistically distributed according to the 
input pulse. Furthermore, this approach allows to elucidate the mechanism leading to energy losses, 
namely a vortical dynamics, which can be advantageously exploited in optimal waveguide design. 
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I. INTRODUCTION 

Optical waveguide design is central to the development 
of efficient devices with applications in optoelectronics 
and nanotechnologies More recently, it has also be- 
come relevant to quantum information processing and 
quantum technologies based on photons [2, 3]. Very of- 
ten these designs involve a complexity that arises from 
the particular details of individual guiding elements as 
well as the existence of a network of waveguides staged 
in series. Thus, to ensure optimal guiding conditions, it 
is of interest to develop efficient and flexible propagation 
and optimization first-principle procedures, more specifi- 
cally, implementing and developing appropriate propaga- 
tion methods to solve Maxwell's equations. These vector 
equations are highly demanding in terms of computa- 
tional effort if adequate balance in accuracy is required. 

For some particular types of waveguides (e.g., optical 
waveguides made of semiconductor materials), though, 
light propagation can be properly accounted for by the 
scalar optics approximation, in particular, the Hclmholtz 
equation [J-|6|- In those cases where the waveguides do 
not bend too severely and the ratio between the inside 
and outside refractive indexes is relatively small, their 
study can be further simplified by appealing to the parax- 
ial or small-angle approximation. The Helmholtz equa- 
tion is then replaced by a simpler scalar wave equation, 
namely the paraxial equation [3, Q . Of course, in more 
adverse cases the paraxial approximation can no longer 
be used and other types of approximations can be con- 
sidered, such as the so-called wide angle equation [7j, but 
this goes beyond the scope of this work. 

The analogy between optical and matter waves is a 
well-known subject in optics [f| lil-fioj. which becomes 
more apparent when considering the isomorphism be- 
tween the paraxial equation and the time-dependent 
Schrodinger equation. Hence, numerical methods avail- 
able in quantum mechanics to propagate wave packets 
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[11| - |13| can also be used to study optimal waveguiding 
conditions for the transmission of light pulses [l4Hl6j - 
actually, in some cases, these methods came originally 
from optics, which shows the fruitful interplay between 
these two branches of physics — , where they are known 
as beam propagation methods (BPMs) . 

Now, apart from having efficient BPMs, one of the key 
questions that arises in a natural way from the quantum- 
optical analogy is whether some extra tools can be de- 
vised and used to analyze and to better understand the 
transport properties inside waveguides. In this regard, 
consider, for example, the analysis of diffraction and in- 
terference with polarized light by means of photon tra- 
jectories EMI, 

which are the direct optical analog of 
the quantum trajectories coming from the Bohmian for- 
mulation of quantum mechanics [20l l2lj]. These studies, 
motivated by single- photon experiments [H|-|24j , allowed 
a better understanding of the flux of electromagnetic en- 
ergy in such contexts, accounting for the formation of 
interference fringes or their disappearance depending on 
the polarization of the interfering beams. In fact, in spite 
of the theoretical nature of these studies, recently pho- 
ton tra ject ories have been inferred from experimental ev- 
idence [251 ] . confirming the topology of those previously 
reported theoretically [l8j |. 

In this work, an analogous concept to that of pho- 
ton trajectory is considered within the context of waveg- 
uiding and utilized to render some light on how trans- 
port is conducted inside waveguides. Unlike the photon 
trajectories mentioned above, here we take advantage of 
the isomorphism between the paraxial equation and the 
time-dependent Schrodinger equation, which allows us to 
readily extend the Bohmian approach to the light trans- 
port accounted for by the former. Since the light flux 
analyzed is time-independent, to avoid possible conno- 
tations associated with the concept of photon trajectory 
(as describing the time-evolution of particles rather than 
fluxes), the concept of optical streamline will be used. 
These streamlines are synthesized "on the fly" in combi- 
nation with the BPM used to obtain the exact evolution 
of the pulse inside the waveguide [lij]. As it is shown, 
the optical streamlines provide us with an alternative, 
complementary insight to the standard wave approach, 
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which allows us to understand the flow of light inside 
waveguides by monitoring its way through locally, i.e., 
just like a tracer particle allows us to determine the flow 
of a classical fluid. In this regard, this constitutes an 
interesting alternative tool to analyze this kind of sys- 
tems and their optical properties. In particular, we show 
how optical streamlines provide us with detailed dynam- 
ical (though stationary) information about guiding prop- 
erties (i.e., the transport of light) and the loss mecha- 
nism. Although there are different designs depending on 
their particular application, in order to properly intro- 
duce the language and methodology, here we focus on the 
Y-junction structure, which is widely used and appears as 
the main component in many other more intricate opti- 
cal devices. More specifically, as a working model, the Y- 
i unction geometry developed by Langer and co-workers 
[26l [27| has been considered, which is closely related to 
realistic waveguides of potential industrial interest. 

This work is organized as follows. In Section [ill a brief 
account on the theory behind both the paraxial equation 
and its connection to Bohmian mechanics is presented. 
Section IIIII deals with the application of this theory to 
the analysis of the light flow through Y-junction- type 
waveguides, providing some numerical results. Finally, in 
Section HV] the main conclusions and future perspectives 
are summarized. 



II. THEORY 



no is the bulk index. Substituting ([2]) into Helmholtz's 
equation ((TJ), the latter can be recast as 
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Here, Vj_ denotes the transverse Laplacian, i.e., the 
Laplacian on the subspace perpendicular to the propa- 
gation direction, = d 2 /dx 2 + d 2 /dy 2 . Apart from 
the paraxial approximation, if we also assume the slowly 
varying envelope approximation holds, i.e., the envelope 
<fi varies slowly in space compared to 2ir/k z , the highest- 
order derivative, d 2 4>/dz 2 , can be neglected in Eq. ([3]), 
which can be expressed as 
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As it can be readily noticed, this equation is isomorphic 
to the time-dependent Schrodinger equation, except for 
the evolution, which is not in time, but in terms of the 
z-coordinate (within this context the space z-coordinate 
acts as the "evolution" parameter). From now on, the 
transversal coordinates (x, y) will be denoted collectively 
by means of the vector R. This quantum-optical analogy 
becomes more apparent if Eq. ([4]) is recast as 
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where 



Consider the optical axis of the waveguide is oriented 
along the z-direction, while its transversal section is par- 
allel to the AY-plane. Moreover, the passage of light 
through the waveguide is describable in terms of a time- 
harmonic electromagnetic field, which allows us to sim- 
plify the treatment by considering a stationary or time- 
independent scalar field $(r) Q. This implies <&(r) sat- 
isfies Helmholtz's equation, 
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plus the corresponding boundary conditions imposed by 
the shape and refractive properties of the waveguide. In 
Eq. n = n(r) is the position-dependent refractive 
index inside the waveguide and k = 2tt/X, with A being 
the light wavelength in vacuum. 

Taking into account this configuration for the waveg- 
uide, in the small-angle or paraxial approximation, the 
scalar field $ can be approximated by a plane wave along 
the z-direction modulated by a certain complex-valued 
amplitude @, i.e., 
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In other words, this means that, for this type of nonuni- 
form waveguides, fast oscillations can be separated from 
slower ones. Accordingly, when moving inside the bulk, 
this splitting allows us to write k z = /3q = kn$, where 



V(v) = ^-[n 2 -n 2 {v)] 



(6) 



is an effective potential function accounting for the 
waveguide refractive profile and 11 = Po = kng plays 
the role of an "optical mass" (2a , [Hj]. Here, in par- 
ticular, we have chosen an optical effective mass to be 
fi = 2/3q in order to enhance the dynamical effects in- 
side the guide. Nevertheless, it is important to stress 
that, while Schrodinger's equation describes a true time- 
evolution, Eq. ([5]) only describes how the stationary, 
three-dimensional scalar field <f> (and, therefore, the light 
intensity) distributes along the waveguide. The fact that 
this field behaves differently along the z-coordinate than 
along the transversal section is used conveniently from 
a computational viewpoint to solve the corresponding 
three-dimensional wave equation, but the field itself is 
always stationary. 

Following with the quantum-optical analogy, consider 
<j) is now expressed in polar form, 
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in order to establish a parallelism with Bohmian me- 
chanics d(J [U , also known as quantum hydrodynamics 
(3fjl l3lj . Within this formulation of quantum mechanics, 
p and S denote, respectively, the probability density and 
phase of the wave function describing a quantum parti- 
cle of mass m, both being real-valued functions. In the 
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present context, though, these quantities are related, re- 
spectively, to the (stationary) energy density inside the 
waveguide and its flux. After substitution of (0 into the 
optical Schrodinger-like equation ([5]) , we obtain a system 
of two real coupled equations, 
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Equation (J5J) is a continuity equation, which describes 
how the energy is distributed throughout the waveg- 
uide. However, given the parametrization in terms of 
z, this equation can also be interpreted as accounting for 
the connection between the energy distribution p = (f>*(f> 
along the transversal section at a certain value of z with 
its flux through this section, J. This flux can be ex- 
pressed as J = pv, i.e., in terms of a local velocity field, 
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Assuming v = EL' = dR/dz and then integrating for x 
and y along the time-like ^-coordinate with different ini- 
tial conditions, one obtains a series of optical streamlines. 
These streamlines allow us to visualize how energy flows 
inside the guide at a local level. That is, we can con- 
nect causally a particular point of the input pulse with 
another one of the output one in a non-ambiguous fash- 
ion — unlike classical rays, the optical streamlines evolve 
in space in accordance with the wave-like behavior of 
the field <f>. In quantum mechanics the corresponding 
quantum streamlines monitor the flow described by the 
quantum (probabilistic) fluid and have been used earlier 
in the literature to describe different physical systems 
[32Tj36j as well as to devise new numerical propagation 
schemes [37| . 

In order to better understand the meaning of two 
points being causally connected, we need to focus on 
Eq. (jSJ). This equation describes the evolution (in terms 
of the z-coordinate) of the phase field S along the waveg- 
uide and presents the form of a quantum Hamilton- Jacobi 
equation 0, [2l|. In it, the effective potential V e g con- 
sists of two contributions, the potential V defined above 
[see Eq. ©] and a kind of "optical" potential, 
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which is analogous to the quantum ■potential of Bohmian 
mechanics. As in the classical version of the Hamilton- 
Jacobi mechanics, solutions of Eq. ^ define trajectories 
as the perpendicular streamlines to surfaces of constant 
phase S, which here translates into the equation of mo- 
tion ([10|), The coupling between © and © through Q 
(or, equivalently, p and its space derivatives) is the reason 
why the wave formulation for waveguides differs from its 
geometric counterpart. The guidance condition (|10[) thus 



allows to properly include typical wave-like features like 
interference and diffraction, for example, in the topology 
of the streamlines, which does not happen, however, in 
the rays one obtains from geometric optics. 



III. NUMERICAL RESULTS 

A. The waveguide design 

Consider a waveguide that transports a light beam in- 
jected at z = to a given distance z = L. Apart from 
its particular geometry, the transport properties of this 
waveguide are characterized by the bulk refractive index 
(no), the guide refractive index (ni), and the cladding 
refractive index (r^)- Physically, these refractive indexes 
can be associated with an effective potential function that 
accounts for the waveguide refractive profile (the profile 
along the transversal direction), as seen in the previous 
section. Inside the waveguide this effective potential or 
refractive profile reads as 



while outside it is given by 
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In order to simplify our study, but without loss of general- 
ity, we are also going to assume that the electric field de- 
pends only on one of the two transversal coordinates, say 
x. This is the condition that determines the confinement 
of light inside the waveguide [6(. The other direction, 
y, will be less relevant, only adding some attenuation 
when moving further away inside the bulk. Thus, in the 
present case, the distribution of light inside the waveg- 
uide is accounted for by a stationary (time-independent) 
scalar field 4>(x, z). Alternatively, taking into account the 
wave picture described in the previous section, (f> can also 
be interpreted as describing the evolution of a pulse of 
light distributed along the x (space-like) direction as the 
value of z (time-like) coordinate increases. 

The configuration of the waveguide along the z- 
direction is as displayed in Fig. [1] It is a Y-junction 
which utilizes a 1x2 multimode interference (MMI) de- 
vice with a lowered refractive index (n w ) in thejunction 
or wedge region, as proposed by Langer et al. [27j (here, 
though, we have assumed the wedge refractive index to 
be equal to the guide one) . This wedge reduces the loses 
over a rather large range of Y-junction angles. Details 
on the numerical values for the geometric and refractive 
index parameters considered here are given in Table [I] 



B. Numerical methodology 

In this work, the Schrodinger-like form of the paraxial 
equation, Eq. ([5]), has been solved by taking advantage 
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FIG. 1: Geometrical layout of the Y-junction design consid- 
ered here, where the significant (geometrical) parameters are 
also displayed (see Tabled]). 



of the computational machinery developed in standard 
quantum wave-packet propagations to solve the time- 
dependent Schrodinger equation. More specifically, we 
have considered the split-operator technique combined 
with the fast Fourier transform method [FIKL3| . As it 
has been shown 0, [IH > this type of BPM results very 
efficient computationally, particularly in optimal control 
scenarios. 

As mentioned above, the purpose of this work is to 
combine the computational performance of BPMs with 
the physical insight provided by the optical stream- 
lines. Therefore, regarding the calculation of the optical 
streamlines, we have chosen to solve Eq. (flQ)) "on the 
fly", i.e., once </> has been evolved from z to z + 5z, the 



TABLE I: Values of the different parameters considered in the 
design of the waveguide used in the calculations shown here 
(see text for details). 



Parameter 


A (fim) 


1 


W (urn) 


2 


Mw {pxn) 


4 


L\ (/nn) 


50 


£ 2 (/mi) 


100 


h dum) 


4 


w (/im) 


2 


Q 


~0.82° 


no 


3.558 


ni 


3.568 


ri2 


3.558 


1 ly_! 


3.568 



streamline is updated, from x(z) to x(z+8z), by means of 
a fourth-order Runge-Kutta integrator fed with the field 
<f>, according to the last equality in (fTQ|) . This is a method 
that has been previously applied to the study of matter- 
wave diffraction by periodic gratings > 38(. Nevertheless, 
from a merely computational point of view, one could also 
consider any of the so-called quantum trajectory meth- 
ods (37|, based on solving also "on the fly" Eqs. ©, © 
and (|10[) . which would skip solving the paraxial equation. 

Once the method is set up, another important element 
when dealing with optical streamlines is that, in order to 
obtain a good representation of the pulse at each value of 
the z-coordinate, one needs to have an optimal sampling 
of the initial pulse, namely p(x, 0). Since one of the pur- 
poses here is to show that the evolution of the pulse along 
the waveguide can be followed by looking at the density of 
streamlines rather than monitoring <f) itself, we have not 
optimized the number of streamlines considered. Thus, a 
total of 11,256 initial positions x$ = x(z = 0) have been 
randomly generated according to the weight p(x,0) [39j . 
This sampling has been made in such a way that the value 
of p(xo, 0) was confined between the maximum of the in- 
tensity distribution (here, at x — 0) and 10~ 3 times this 
value, i.e., p(x — 0, 0) > p(x , 0) > 1CT 3 x p(x = 0, 0), for 
all xq [see Fig. Ufa)]. These initial conditions have then 
been propagated and, at each value of z, a histogram 
was built up with them to compare with the intensity 
obtained from the paraxial equation. This step in the cal- 
culation is crucial in order to obtain a fair reproduction 
of the optical result. We have observed that the calcula- 
tion of the streamlines "on the fly" , as described above, 
does not decrease much the efficiency of the BPM consid- 
ered when computing a low number of sampler stream- 
lines, i.e., streamlines calculated to determine the flux of 
light without caring much about having a good sampling. 
Otherwise, the time of computation increases strictly lin- 
early with the number of streamlines, as in any standard 
trajectory calculation. 



C. Analysis of the results 

In Fig. [UJa), it is observed how the histogram fairly 
reproduces the profile (along the waveguide transversal 
section) of the input pulse at z = /im. The same is 
also found for larger values of the z-coordinate, at z = 
450 /im, as seen in Fig. [2jb). In general, this behavior 
can be found at any value of z, which allows us to study 
the transmission properties along the waveguide by only 
considering the propagation of the histogram, as shown 
in Fig. [3l In this plot, the intensity reaches a maximum 
as soon as it leaves the 1x2 MMI device and then splits 
up into two identical fluxes, each one displaying a series 
of bounces inside its corresponding waveguide as they 
propagate. These bounces can be somehow related to 
the bounces that a light ray would undergo according 
to a geometrical optical viewpoint (the ray propagates 
along a straight line until it reaches the boundary of the 
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x (urn) x (urn) 

FIG. 2: Initial (a) and final (b) intensity transmitted through the waveguide. In each panel, black squares denote the optical 
streamline calculation (histogram), while the red solid accounts for the value obtained from the standard optical calculation 
(paraxial equation). The parameters used in these calculations are given in Tabled] 



waveguide, which defines its turning point in accordance 
with the laws of reflection and refraction) . 

The analogy with geometrical rays can be extended 
by means of the optical streamlines here defined. In 
Fig. 21 a set of 200 of such streamlines, chosen (in pro- 
portion) from among the full set run, are plotted. As it 
can be noticed, optical streamlines behave as geometrical 
rays, undergoing a series of bounces and reflection at the 
boundaries of the waveguide, although their flow is more 
laminar, like tracer particles moving on a fluid. This is 
a consequence of the fact that they follow a wave and 
therefore they have to accommodate to the features un- 
dergone by such a wave, such as diffraction, interference 
or refraction. In connection with this property, also note 
that they satisfy the non-crossing property of Bohmian 
trajectories fl(j : two streamlines cannot pass through the 
same position x for the same value of the evolution co- 
ordinate, z. This allows us to establish a clear difference 
between streamlines evolving along the two branches of 
the Y-junction. In the case considered here, this might 
seem to be not so important, for the evolution along z 
is symmetric with respect to x = 0. However, notice 
that as soon as a slight asymmetry between pathways 
would appear (different refractive index, length, shape, 
etc.), the splitting would also be different. In such a 
case, the optical streamlines will provide us with unam- 
biguous information on which specific portion of the in- 
put signal goes through each branch. Thus, this type 
of information, which cannot be disentangled in wave- 
based treatments by their intrinsically global nature, can 
be of interest in applications based on waveguide tech- 
nologies, such as Mach-Zehnder interferometry or other 
optical interference devices. Actually, even in the present 
case, this propert y re sults useful, because it allows us to 
causally connect [4l| a certain feature from the output 
signal at z — 450 /im with a particular region of the in- 
put, at z = 0. Later on this information can be used 
to make the sampling only over a certain subset of ini- 
tial conditions (e.g., those that does not lead to losses at 
z = 450 /urn). 



Apart from the properties above, an additional in- 
teresting feature enabled by an analysis based on opti- 
cal streamlines is the detection of energy losses through 
the waveguide boundaries. Detection of such losses by 
standard means, i.e., by only studying the evolution of 
<p(x, z), the intensity p(x,z) needs to be integrated over 
the waveguide cross-section (at a given z value). From 
a more graphical viewpoint, the intensity should be en- 
larged until the amount of intensity that is flowing out- 
side is perceived [see Fig.@Jb)]. For example, in Fig.[^b) 




FIG. 3: Three-dimensional representation showing the evo- 
lution of the initial pulse p(x, 0) as it propagates (along the 
z-coordinate) throughout the waveguide (below, the corre- 
sponding contour-plot is also displayed). In particular, p(x, z) 
has been obtained as a histogram-like distribution of optical 
streamlines at each z value, which corresponds fairly well with 
the representation one would find by standard wave propaga- 
tion methods (see text for details). The parameters used in 
these calculations are given in Table [I] 
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FIG. 4: Optical streamlines (a) and contour-plot of the evolution of the histogram (b) illustrating the flow of light throughout 
the waveguide. In both figures, the straight black lines mark the waveguide boundaries (see Fig. fT}. In part (a), blue and red 
lines are used to emphasize the fact optical streamlines also satisfy the non-crossing rule that characterizes both quantum [iO ] 
and optical fluxes [18J. In part (b), the maximum value of the contours is set up to 0.1 in order to appreciate the outgoing light 
flow. White labels correspond to different vortices that appear during the evolution of the flow (see text for details). 



the presence of two "wings" lying out of the waveguide 
boundaries is clearly seen, which indicate the leakage 
of light from the waveguide. However, they are neither 
perceived in the contour-plot of Fig. [3J nor the continu- 
ous outwards flow of energy (or how and when this flow 
happens). A more convenient and insightful way to un- 
derstand this flow arises by observing the topology of 
the optical streamlines, as it is apparent from Fig. EJa). 
Although the outermost streamlines are associated with 
very low intensity (that is precisely the reason why they 
appear so sparse in the figure), it is readily noticed how 
the flow is outwards. Actually, whenever a node of the 
intensity appears, it pushes away the streamlines, as can 
be seen at z « 150 /im, 275 /im or 450 /im (and inwardly 
at z rs 400 /im). In this sense, optical streamlines ac- 
quire the role of optical tracers of the energy flow inside 
the waveguide, allowing us to monitor it at each step. 
This property thus results very useful in the optimal de- 
sign of waveguides and the control of light inside them, 
for it gives us a straightforward method to test losses and 
gains. 

At a more detailed level of description, the optical 
streamlines provide us with a better understanding on 
the dynamics that is taking place inside the waveguide, 
such as the different features which one can observe later 
on in the modes or energy losses. In this regard, consider 
first the high intensity peak observed at the entrance of 
the second section of the 1x2 MMI (see Fig. [3j, which 
later on gives rise to a node or vortex [V w , in Fig. HJb)] 
in the wedge center. By inspecting the optical stream- 
lines [see Fig.@|a)], one readily notices how most of them 
quickly undergo an inwards motion as z increases, just 
before V w becomes incipient at the position where the 
second section of the 1x2 MMI starts (regarding the 
streamlines that surmount the vortex through the op- 
posite side; see comments below). By construction, the 
tight squeezing of the swarm of streamlines takes place 
just at the entrance of the second section of the 1x2 
MMI, from which it undergoes a fast expansion. This 



expansion causes the splitting of the beam into two dis- 
tinctive halves (partial beams), each one channeled into 
one of the branch of the Y-junction and, moreover, a vor- 
tex in between. It is interesting to compare this behavior 
to that of a classical fluid which is released from a pipe, 
thus giving rise to channeled streams and hollows or low 
density regions in between. 

Now, consider only one of the branches of the waveg- 
uide, for example, the left-hand side one. From a simple 
inspection of Fig. @|b) , where the contours refer to low 
intensity values (up to 10% of its maximum, at the en- 
trance of second section of the 1x2 MMI; see Fig. [3]), one 
observes that the flow of energy appears and disappears 
alternatively along two pathways as z increases (see black 
dashed lines). This is an effect of the modes generated 
within the branches of the waveguide that can be easily 
explained in terms of a simple model. Consider that two 
kind of wave packets travel along each pathway. They 
are being populated alternatively as they evolve along z, 
i.e., when one starts to gain population, the other losses 
it, just as in a simple two-level system, 

4>{x, z) = cos(uJz)tfi(x) + sm(ujz)(p2(x)- (14) 

According to this simple model one can easily explain 
that, for example, near the vortex Vi [see Fig. BJb)], one 
of the pathways starts to get depopulated very quickly, 
while the other reaches its maximum (see also Fig. [3J. 
The same behavior is again observed when approaching 
V3. However, although this two-level model describes the 
features in the population changes as the optical modes 
travel through the waveguide branches, no information 
is provided about the population transfer between path- 
ways or energy losses. These are processes suitable to be 
explained by means of the associated optical streamlines, 
for they provide us with precise information about the en- 
ergy flux and its evolution along the waveguide, as seen 
above. Thus, as it can be seen as the mode approaches 
V2 [see Fig.QJa)], the streamlines start to move smoothly 
rightwards, from one pathway to the other, until their 
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number has decreased considerably in the first pathway. 
Near V 2 there are three options for those streamlines still 
remaining along the first pathway: (1) some of them un- 
dergo a sudden turn and will enter into the second path- 
way stream, (2) some will slightly bend and will con- 
tinue along the first pathway beyond the vortex, and 
(3) another group leave the waveguide after the bend- 
ing, thus giving rise to the losses (eventually, streamlines 
from group 2 will also scape). This behavior can be ob- 
served not only around , but around any other external 
vortex (Vi, V 3 and the homologous vortices along the sec- 
ond pathway), i.e., any vortex facing the boundaries of 
the waveguide. The goal of waveguide optimal design is 
to avoid losses and therefore it can highly benefit from 
this optical streamline representations. 

IV. SUMMARY AND FUTURE PERSPECTIVES 

The purpose here was to introduce the concept of opti- 
cal streamline in the literature of waveguides as an ana- 
lytical tool, in combination with efficient BPMs to carry 
out their propagation. To carry out our hydrodynam- 
ical analysis, we have considered as a working model 
a Y-j unction with a geometry closely related to realis- 
tic waveguides of potential industrial interest, originally 
proposed some year ago within the broad family of Y- 
junctions [27| . The particularity of these devices is that 
their design seeks for a minimum loss of the input while 
having a clean splitting of the signal, at the same time. 
Thus, given the interest and relevance of these designs in 



modern technologies, to have at our disposal additional 
tools which can help us to design, study, explain and un- 
derstand such systems is of great importance. In this 
regard, here we have mainly focused on determining how 
the modes travel through the wave as well as on eluci- 
dating the energy/light loss mechanism in waveguides. 
These two features can be efficiently controlled by means 
of an appropriate design [15l . l26j , which is precisely the 
scope of the waveguide optimal design techniques and 
the eventual target of our optical streamline analysis, 
where it can be of much interest and help. This is ac- 
tually a natural and important extension of the present 
work, which is currently under development. Nonethe- 
less, we also like to stress the general usefulness of the 
proposed new tool, which might be of relevance in other 
more complex scenarios where interference or merging of 
different signals could be more efficiently studied (e.g., 
Mach-Zehnder interferometers [IB])) as wen as other do- 
mains going beyond the paraxial approximation (e.g., the 
wide-angle equation [7J). 
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